Effects of gating and contact geometry on current through conjugated molecules 

covalently bonded to electrodes 
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We study the effects of gating and contact geometry on current through self-assembled monolayers 
of conjugated molecules strongly coupled to gold electrodes by sulfur "anchor groups". The current 
changes by more than an order of magnitude depending on the angle between the axis of the benzene- 
dithiolate molecules and the normal to the electrode on the less coordinated "top site" position. 
The effect of gating is also much stronger in this case compared to higher coordinated "hollow 
site" binding of the molecule on a Au(lll) surface. The large hybridization of the molecular states 
with electrode states for the hollow site leads to practically ohmic current-voltage characteristics. 
Changes in molecule-electrode geometry accompanying the gating of the SAM may be the reason 
for strong changes of the conductance observed by Schon et al. in the "slot" geometry. 

PACS numbers: 73.40. Gk, 73.61.Ph, 73.63.Rt 



I. INTRODUCTION 

Studies of electron transport 

through organic molecules (molecular films), viewed as 
the possible components of molecular electronic devices, 
are a very active area of research Although the 

rectifying properties of the molecules in two-terminal de- 
vices were demonstrated in 1990s [Q, only recently have 
f/iree-terminal devices based on a self- assembled mono- 
layer (SAM) of simple conjugated molecules containing 
either one or two phenyl rings CgHg or two or three thio- 
phen rings C4H4S been fabricated The molecules 
were gated from the edge of the vertical device through 
L =30 nm of Si02 gate oxide (30 nm was the nom- 
inal oxide thickness, but the high-resolution transmis- 
sion electron microscopy showed much thinner oxide bar- 
rier of 4-5nm ||^). In the slot geometry of this exper- 
iment, the drain current Id through only several thou- 
sand molecules close to the gate oxide should have been 
affected by the gate voltage Vg. The measured con- 
ductance dId/dV at zero source-drain bias was about 5 
/iS/molecule. The gating field on the molecules should 
be small because the length of the molecules (the "chan- 
nel" length) is tiny, 2t = 1 — 2 nm, so that the geo- 
metrical aperture factor is small, t/L <^ 1. In spite of 
this, Schon et al. have initially estimated a dramatic 
change of conductance of the molecules by a factor to be 
about 10^ (it was later determined that the effect was 
overestimated, probably because of disorder and other 
factors in pure SAMs |^]), with a large transconductance 
dId/dVg = 12-13 mA/V for pure SAMs of biphenyl 
molecules and about 6 mA/V for phenyls [|j at room 
temperature. Dilute systems of conducting molecules 
(namely, 4,4'-biphenyldithiol and some other species) em- 
bedded in a matrix of non-conducting alkanedithiols have 



also been measured with the aim to sample individual 
conducting molecules Schon et al. have found that 
the conductance through the molecule is about 0.1/iA/V 
at zero drain voltage V = 0, and approaches one con- 
ductance quantum 2e^//i — 77.5/iA/V at relatively small 
gate voltage Vg « —0.3 V. The histograms of conduc- 
tance measured on a series of devices suggested that the 
conductance through the gated molecules was quantized 
in the fundamental units of 2e^ /h. The conductance peak 
position varied with the gate voltage Vg from device to 
device. The transconductance dId/dVg was found to be 
in the range of 150/iA/V at 4K, decreasing to about 
IQpiA/V at room temperature (in the earlier measure- 
ments the peak value of dId/dVg was estimated of about 
13mA/V 1^ for the biphenyl molecules). The correspond- 
ing modulation of conductance in the film of individual 
dithiol molecules has been estimated to be much larger 
than 10"^. This is a huge effect indeed if it is to be ex- 
plained by a mechanism related to the behavior of indi- 
vidual molecules. 

Very recently this work has been extended to SAMs of 
benzene-(l,4)-dithiolate (-S-C6II4-S-), the simplest con- 
jugated molecules, referred to below as BDT, with only 
one phenyl ring . The area of the measured device com- 
prised about 10^ molecules, mostly non-conducting alka- 
nes with a small fraction (10""*) of the conducting BDT 
molecules 0. In this geometry, only one or a few BDT 
molecules should have contributed to transport. The 
first peak in BDT conductance has been found at a bias 
of about 2.1V with the corresponding value of 25/iA/V 
0. Interestingly, the conductance appeared to be much 
larger than that found in the earlier break-junction ex- 
periments by Reed et al. ||] who observed the first peak 
in conductance of 0.05/iA/V at 1.4V. Schon and Bao have 
also observed a jump in conductance from 0.7-1.0 /iA/V 
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to much larger values close to the conductance quantum 
77.5 /xA/V at the gate voltage Vg w —0.4 V and the drain 
voltage of about —0.7 V (this is a jump of almost two or- 
ders of magnitude) . The position of this switching shows 
some hysteresis when the gate voltage is varied in the 
range —0.6 < Vg < V. This hysteretic behavior can- 
not possibly come from the individual molecule behavior 
and should be related to some extrinsic mechanism, as 
we shall discuss below. Also interesting is the observed 
change in the fine structure of the conductance: there is 
a series of small peaks 45 meV apart from each other on 
the conductance curve for = 0, and the period shrinks 
down to 40 meV at = —0.5 V. Such peaks are usually 
attributed to tunneling assited by molecular vibronic ex- 
citations, so the change in the period indicates a possible 
change in molecular conformation. These abrupt changes 
cannot be explained by the electric field created by the 
gate, since the corresponding field on the molecule should 
be small. 

With regards to a possible origin of the observed be- 
havior, we first mention the importance of the geometry 
of the molecule-electrode contact |^ . If the orientation 
of the molecule with respect to an electrode changes, so 
does its conductance. Indeed, we have found earlier that 
the conductance of BDT (or any other conjugated thiol- 
terminated molecule) strongly depends on the angle 9 
between the molecular "backbone" and the normal to 
the gold surface (Fig. |l|). In a simple "toy" model, the 
current dependence is oc sin^ 9 in the regime of strong res- 
onant tunneling (large bias) and even stronger, c>c sin"* 9 
in the regime of non-resonant tunneling (small bias) . 
Note that the angular dependence of current is much 
stronger when the end sulfur is in the less-coordinated 
"top-site" position above a surface Au atom, compared 
to a "hollow-site" position, when it is bonded to three 
surface gold atoms, see Fig. |^. These bonding positions 
were considered in the literature as being the most fa- 
vorable Indeed, it has been provisionally suggested 
that the conformational changes of the BDT with respect 
to electrodes, discussed in are responsible for the ob- 
served "switching" behavior of SAMFETs 0. Also, one 
cannot exclude that the observed changes in conductance 
through BDT are caused by charge trapping-detrapping 
processes close to the interface between the molecular 
layer and the gate Si02. 



II. SELF-CONSISTENT CALCULATIONS OF 
ELECTRONIC STRUCTURE AND TRANSPORT 

In order to gain more insight into the origin of the 
observed extraordinary dependence of molecular conduc- 
tance on a weak gating field observed in Refs. Q,^, we 
have performed a series of self-consistent calculations for 
different attachments of the BDT molecule to Au(lll) 
surface. We have found that the effects of charge redis- 
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FIG. 1. Schematic of the benzene-dithiolate molecule on 
top and hollow sites. End sulfur atoms are bonded to one and 
three surface gold atoms, respectively. 9 is the tilting angle. 



tribution as a function of molecular configuration and/or 
external field are important, since the mismatch between 
the work functions of Au and BDT results in a consid- 
erable charge fiow between the molecule and gold elec- 
trodes. Without external bias, the lowest unoccupied 
molecular orbital (LUMO) lies closer to the Fermi level 
of Au compared to the highest occupied molecular orbital 
(HOMO). The current through the molecule strongly de- 
pends on both the tilting angle and the self-consistent 
charge redistribution across the molecule. Although it 
is not clear a priori what the sign of the combined ef- 
fect would be, we have found that our previous conclu- 
sion about strong orientation dependence of the conduc- 
tance through anisotropic conjugated tt- orbitals and the 
s-orbitals on electrode Au remains valid 

BDT attaches strongly to the gold substrate by thi- 
olate end groups -S- that form covalent bonds with Au 
|p|,pl| . In order to properly account for such a bonding in 
the present calculations, the Au atom(s) connected to S 
are treated separately from other gold atoms that com- 
pose the gold electrode. The conductance is computed 
with the use of the general procedure of Ref . Q . The 
gold electrodes are described by a tight binding model 
with nine s— , p~, and d— orbitals per each Au atom 
with parameters from Ref. p^ ] . The equilibrium molecu- 
lar geometry is found by total-energy density functional 
minimization p^ . The tight-binding parameters for the 
molecules and molecule-lead interfaces are taken from the 



solid-state table of elements [IS 



The onsite energies in the present tight-binding model, 
which are very important for finding the correct charge 
redistribution between the molecule and the electrodes, 
have been estimated from the Hubbard model in the 
atomic limit. The energy of an isolated atom is approx- 
imated as Em = Eq— eAqm + ^UAq^. Here Eq is the 
energy of a neutral atom with the atomic energy level at 
— e < with respect to the vacuum level (energy origin), 
Aq is the excess charge of an atom, and U is the in- 
traatomic Coulomb repulsion. In this approximation, we 
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obtain e = ^ {A + I) , U = I — A, where A and / are the 
experimental atomic values for the affinity and the ion- 
ization energy, respectively. These expressions have been 
used to estimate the following parameters used in the 
present work: U — 11.5, e = 7.8 for H; U — 6.3, e — 5.2 
for C;U = 7.8, e = 6.5 for S; and U = 6.7, e = 5.9 for 
Au (all in units of eV). We would like to mention that 
the use of different values for the one-electron energies e, 
like the ones from Ref. [^) leads to unphysically large 
charge transfers. 

We have calculated the current through BDT on 
Au(lll) in both the top-site and the hollow-site posi- 
tions. Including onsite and intersite Coulomb interac- 
tions one finds that the onsite one-electron energies for 
state a at the site m should be adjusted as 



(1) 



where ej^^j are the onsite energies in a system with neutral 
atoms, Agm are the charges on sites, 7mm' = l/|m — m'|, 
(/)^ the image potential, and e < is the electron 
charge ||l^. The charge Agm is found self-consistently 
from the local density of states, which is given by the 
site-projected imaginary part of the exact Green's func- 
tion of the problem. The total retarded Green's func- 
tion Gmam'a'{E) is Calculated by "attaching" the semi- 
infinite leads to the molecule |l^. As a result of the 
attachment the molecular levels acquire a width that 
strongly depends on the coupling between electrode and 
the molecule, F ~ ^au-s/^Au- Here Dau is the width of 
the s-band for Au electrodes, ^au-s is determined mainly 
by the spa hopping integral from Au to the end sulfur 
atom on the BDT molecule, which is of the order of 1- 
2 eV. One should expect significant broadening of the 
molecular levels when the molecule is attached by a thiol 
group to Au, since a strong chemical bond is formed. 

Under zero bias voltage, the electron charge on the 
site m can be found from the Green's function in the 
standard manner as 

/oo 
dEN^a{E)f{E), (2) 

iVma (E) = ~ -ImGraama (E), (3) 

where N,na{E) is the density of states a on the site m, 
f{E) = (l -I- exp ^Zj^" ) ^ is the Fermi function, and Ep 
is the chemical potential found from the global charge 
neutrality of the system. 

In the case of finite bias voltage the system is out 
of equilibrium and one has to find the charge that is 
"flowing-in" from the electrodes onto the molecule, cf. 
p6[ . Then the DOS on the site m, related to the influx of 
electrons from the lead w = 1, 2..., is written as N^{E) — 



2(for spin)X;fe™k 
where V'm {kz^^^-w 



(fc!f,k|]')|%(i?-£;^^k„) 



is the wave function at the molecular 
site m, which asymptotically becomes an incident Bloch 
wave in the lead w far from the molecule with the wave 
vectors (/c2,k||). Now we can find the occupation num- 
ber for that site on the molecule due to charge flowing 
in from the lead w as = E» IZ. ^EN^ (E) U{E), 
where fw{E) is the Fermi function for the wth. lead (i.e. 
with Ep = Ep^). In order to calculate the Green's func- 
tion (and the charges Qm) we deflne the "channels" such 
that kz — k'^i{E), where / = 1,M enumerates all the 
quantum states in the lead unit cell (slice) ||l^. It is 
convenient to re-write the expression for the charges in 
terms of "open channels" . The "open channel" is defined 
as a Bloch wave which propagates in the lead at a given 
energy. The Bloch waves incident on the molecule (i.e. 
having the velocities towards the scatterer, vi > 0) wfll 
contribute to the charge flowing to the molecule from a 
particular lead: 

NZiE) = 2j2l- dkl I V^m (A^r , k|[) ^5(e- El 



(4) 



where dz is the unit cell length along the lead. V'm are 
normalized for the length of the wire, which drops out of 
the final answers. Note that the integration in (^ goes 
over the whole Brillouin zone, not just over kz > 0. The 
delta- function picks up the open channels on the leads. 
From now on we can drop the lead index and assume that 
one can later sum up all the charges flowing from all the 
leads. Once the Hamiltonian is set up, one calculates the 
charges on the sites, recalculates the onsite energies Cma 
and continues iteratively until the charges converge. 

The current through the film is given by a standard 
expression ]l7| , p^ 



dE 



T(i?), (5) 



where q = \e\ is the elementary charge, and T{E) is the 
transmission probability 



T{E)= Vnn'{E,)<i 



(6) 



where the summation goes over the surface Brillouin zone 
of the lead. Transmission coefficients i„„' {E, k|| ) between 
the scattering channels n and n' are found from the so- 
lution of the scattering problem |jl^ . In the case of weak 
molecule-electrode bonding the transmission probability 
is approximately given by the Breit-Wigner formula |9| 



T{E) 



E 



TrjX 



rii rR 



[E - ErY + {TrL + r.fl.) 2/4 



(7) 
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where enumerates the energies of the molecular or- 
bitals (MOs) contributing to transport (not all of them 
do, see Figs. ||, TrLiR.)/^ is the rate of the carrier 
transfer to the left (right) electrode from the molecu- 
lar orbital r. This formula applies when the width of 
the MOs is much smaller than the energy difference be- 
tween them, so that the resonances do not overlap. Each 
conducting molecular orbital produces a step-like contri- 
bution to the current. Indeed, when the resonance falls 
into the "window" between the lowest and the highest 
Fermi levels in the leads EpL < Er < EpR, the current 
obtained from Eq. || is 

~ TF — ■ 

n i rL + i rfl 

It follows from this analysis that the current-voltage 
characteristic should look as a series of steps, occurring 
when the resonant conditions are satisfied for particu- 
lar conducting molecular orbital. The apparent negative 
differential resistance (NDR) at bias above 2V, Figs. ^j6[ 
results not from resonant tunneling but from the elec- 
trode density of states. In the present model the elec- 
trode DOS in bounded from above for each particular 
value of k|| . As a result, the current will be zero at the en- 
ergies above some threshold. The apparent NDR persists 
in the present calculations irrespectively of the number 
of basis functions (s-, sp-, or spd-basis), Fig. ^. 

III. EFFECTS OF CONTACT GEOMETRY AND 
GATING ON CURRENT- VOLTAGE 
CHARACTERISTICS 



The s-electron on the top Au can only hop onto a S z- 
orbital via a {z\H\z) = ssa hopping integral. Obviously, 
for BDT on the top site, this result holds for all incident 
electrons with any k|| . Therefore, the x and y p-orbitals 
on S cannot be traversed by electrons incident from the 
contact. At the same time, only those states on the sulfur 
ion are coupled to conjugated tt— orbitals on the benzene 
ring. Therefore, for the BDT on the top site and ori- 
ented normal to contact, the current will be suppressed, 
as observed in calculations by Di Ventra et al. p8| , p^ . 
Obviously, this symmetry selection rule is lifted for any 
9 ^ Thus, the previous calculations ||l^,|l^ have been 
performed at an artificial singular point. Incidentally, 
the same conclusion applies to the scattering of the car- 
riers incident with k|| ~ (surface F— point) on upright 
BDT on a hollow site. Indeed, in this case the matrix 
element for hopping to the sulfur atom on the molecule 
is proportional to '^i{xi\H\s) exp(ik||/9i) oc = for 
k|| = 0, where U are the directional cosines connecting 
the center of the triangle formed by three Au atoms on 
Au(lll) surface with Au atoms in the corners at posi- 
tions Pi. The same is obviously true of the hopping to 
the y p-orbital on S. Thus, in the case of BDT on a hol- 
low site, all the contribution to the current comes from 
states with k|| ^ 0. Therefore, for BDT placed upright 
on the hollow site, the total current is not suppressed, as 
it is for BDT on the top site. Consequently, the current 
for BDT on the hollow site is considerably less sensitive 
to the precise contact geometry. 

A. Density of states and transmission 



We argue below that gating of SAMs in the exper- 
iments [|7y2l[ may have led to changes in molecule- 
electrode geometry and, consequently, to large changes 
in conductance. Given the strong orientational depen- 
dence of the current through conjugated molecules like 
BDT, and that in experimental SAMs the molecules are 
never positioned strictly normal to the electrode surface 
(as was assumed in Refs. p^ , p^ ), we shall present the 
results for the transmission, density of states (Figs. ||, H) 
and I-V curves (Figs. for a series of tilting angles 
9 between the backbone of the molecule and the normal 
to the Au(lll) surface. The 0— dependence of the I-V 
curves for BDT on the top site and hollow site is illus- 
trated on Figs. ^,|| for 6* = — 30°. It is especially strong 
for BDT on the top site. The majority of the results 
is given for 9 — 10°, which seems to be a reasonable 
choice for experimental SAMs. Note that in the upright 
position = in top site (i.e. perpendicular to the con- 
tact surface, as was assumed in Ref. [|l8|) the overlap 
between the S x and y p-orbitals {xy being in contact 
surface plane, z normal to the contact) and the s-orbital 
on Au (or jellium) is exactly zero by symmetry, since 
{x\H\s) = {y\H\s) = 0, where H is the Hamiltonian. 



Most of the present results can be appreciated from 
the analysis of the transmission probability T{E) and 
the density of states N{E) on the BDT molecule, see 
Figs. ||,||. One expects from the golden rule that the 
transmission would be proportional to the density of 
states. However, although peaks of both functions fol- 
low each other rather closely, there are important differ- 
ences between the density of states and the transmission. 
It is easier to analyze the results for the top position 
first. There are two sharp peaks around the Fermi level 
Ep, marked as tt* (at Et^, = Ep + 0.5 eV) and tt (at 
-Ejr — Ep — 1.0 eV), Fig. ^(a). Transmission is almost 
zero at i5 > Et^* , but there is a large density of non- 
conducting states in this energy interval. Those non- 
conducting states are formed at the end of the molecule 
and reside primarily on the end sulfur atom and gold 
atom on top of which the molecule sits, with little cou- 
pling to Cvr conducting states on the ring. The tt* peak 
contains mostly Au and S states and some Ctt ring states 
whereas the tt peak is made mostly of S and Ctt ring 
states with a little addition of Au states. Sulfur atoms 
introduce the states in the HOMO-LUMO gap of the 
benzene ring (which is about 6.5 eV) and hybridize with 
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FIG. 2. Density of states N{E) and transmission T{E) 
through benzene-dithiolate (BDT) molecule on Au(lll) as a 
function of energy: (a) BDT on the hollow site, (b) BDT on 
the top site (see text for the description of the configuration). 
The broken line indicates the transmission T{E) under the 
bias voltage 2V. Molecules in both configuration are tilted by 
10°. 

Ctt states to make the conducting pathways across the 
BDT molecule. As a result, we see a much smaller gap 
between tt and tt* states in BDT, which is iJjr* — -E-tt = 1.2 
eV. In the case of the hollow-site position, the situation 
is considerably different. There the tt and tt* states are 
much broader and pushed apart by much stronger hy- 
bridization with three underlying Au atoms than for the 
top site. Fig. §(b). The "soft" energy gap for the hollow 
position is E^, -E^= 3.25 eV, Fig. | (a) (cf. 

It is instructive to compare the present results with the 
recent data for Ni-BDT-Ni SAMs [|l| and jellium-LDA 
calculations jist . The systems should have many similari- 
ties since Au and Ni have almost the same work function 
of about 5.1 eV In both systems a large conduc- 

tance peak is observed at a bias of about V = 2.1V, 
which is larger than the earlier value of 1.4V The 
peak in Ni system is narrower and does not show any ap- 
preciable spin splitting. Interestingly, in Ni-BDT-Ni the 
additional resonant-like features are found at low bias 
voltages of about 0.3 V and 0.9 V. They correspond to 
smaller conductance compared to the peak at 2.1V. The 
position of the peaks should be compared with that of Ni 
d-states. It is well known that the energy of a minority 
peak in Ni DOS is very close to the Fermi level, whereas 
the majority peak is about 0.5eV below Ep- If the con- 
ducting molecular states were considerably smeared out, 
like in the case of hollow position. Fig. g, then the ob- 
served low-bias features might be due to those peaks in Ni 
d-DOS. The fact that the first peak is observed at 0.3V, 
and not much smaller bias, may be due to energy de- 
pendence of the molecular density of "tail" states at Ep, 



which shifts the peaks by 0.3V. If, on the other hand, 
the peaks in the molecular density of states are sharp, 
like in the case of top position. Fig. |2[ then the peaks 
in conductance should correlate with the position of the 
molecular orbital closest to the Fermi level (LUMO, ac- 
cording to the present work). The position of the LUMO 
in BDT in the present model is about 0.5 eV above the 
Fermi level for the top position. Therefore, we expect 
that in the top-site configuration there should be two 
peaks in the conductance, one at about 0.5-1.0 V (the 
position of the d-peak in Ni DOS with respect to the 
TT* resonance, depending on the connection between the 
molecule and the electrodes) and another at about 0.5 V 
higher than the first one (at 1.0-1.5 V). Interestingly, this 
is very similar to what was observed experimentally for 
the Ni-BDT-Ni system, with the peaks &tV = 0.3 V and 
0.9 V [in. If, on the other hand, the LUMO (tt* state) is 
at leV above the Fermi level, than the spin-peaks in Ni 
d-DOS might have produced the peaks in conductance 
at 2-3V. Indeed, there is a conductance peak at 2.1V in 
both Ni- and Au-based systems, but it is not spin-split 
in the case of Ni electrodes. It is worth mentioning that 
the position of the first peak in conductance is a strong 
function of the tilt angle, and it may vary significantly, 
Fig. |. Finally, we note that the HOMO-LUMO gap in 
the "jellium" calculations is 5 cV [|l^, which is sub- 
stantially larger than that in bare BDT molecules, and 
that is unlikely. The calculated value of the first peak 
in conductance in jellium-LDA is 2.4V, larger than the 
observed value of 2.1V. This suggests that the calculated 
gap is larger than the measured one, which is contrary to 
the usual notion that the LDA gaps are smaller than the 
experimental ones due to insufficient account of electron- 
electron correlations. As follows from this discussion, one 
needs more analysis to draw definitive conclusions about 
the position of the lowest conducting orbital with respect 
to the Fermi level in electrodes in Au-BDT-Au and Ni- 
BDT-Ni. 

B. Gating the molecules 

The gating effect on the transmission and I-V charac- 
teristics is illustrated in Figs. The gating is modeled 
by shifting the on-site energies on the molecule by ^g, 
which is usually —0.5, 0, and 0.5 eV in the calculations. 
Obviously, in the experimental situation Q such a large 
shift would require very large gating fields, comparable 
to the atomic fields in the order of magnitude. This is 
because one has to substantially change the electronic 
states on the molecule, and the characteristic energy is 
given by the HOMO-LUMO gap, usually a few electron- 
volts. Such large fields could not be possibly produced in 
the slot geometry with the channel length of only i = 1 — 2 
nm through the gate oxide with thickness L = 4 — 5 nm 
0. Schematic of this gate is shown in inset in Fig. ||(a). 
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FIG. 3. Transmission T{E) through the BDT molecule 
on Au(lll) and the effect of gating: (a) hollow site, (b) top 
site. The gating is simulated by shifting the on-site energies 
on BDT by $3 =0.5, 0, and -0.5eV. Note the presence of a 
sharp peak in T{E) for top site (originating from the LUMO 
on BDT) close to the Fermi level Ef- 
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FIG. 4. Current-voltage characteristic and effects of gat- 
ing and tilting with respect to the Au(lll) electrode surface 
on current through the BDT molecule on the top site: (a) 
effect of the gating, onsite energies are shifted by the amount 
"l>g indicated on the figure, (b) effect of increasing tilt angle 
9. Current is in units of Jo = 77.5fj,A. Inset: schematic of 
the gate geometry, 2t —1 nm, L = 4 — 5 nm (nominal oxide 
thickness was 30nm), Vg is the gate voltage. 



The analytical solution to this electrostatic problem can 
be found by standard methods and it naturally contains a 
small parameter t/L ^ 1, so the gating on the molecule 
itself would be much smaller than the nominal gating 
voltage Vg. 

One can speculate that large gating may result from 
e.g. charge accumulation in the gate oxide next to the 
molecular film. However, changing the oxide from Si02 
to AI2O3 apparently has not modified the results 0. Be- 
sides, there is an abrupt change of conductance by about 
an order of magnitude at the gate voltage Vg — —0.3V, 
which would suggest a high sensitivity of the interface 
charge to the bias voltage. Both facts are difficult to rec- 
oncile with the idea of interface charge accumulation but 
we will study this possibility. 

Comparing Figs, ^(a) and ||(b), one observes that the 
smaller hybridization between S p- and Au s-states for 
BDT on the top site produces sharp features in the energy 
dependence of the transmission at the Fermi level. The 
LUMO in this case is above the Fermi level by only about 
0.5eV. The shifts of on-site energy by similar amount 
substantially change the transmission at the Fermi level. 
Fig. ^(a), and the corresponding current per molecule as 
shown in Fig. |^(a). There is a pseudogap at low volt- 
ages V ^ IV, with the threshold voltage moving by an 
amount comparable to the external shift <i>g for the top- 
site configuration. By contrast, the large hybridization 
of S p-states with Au on the hollow site results in much 
broader energy tails of the resonant peaks in the density 
of states in the gap region in the vicinity of the Fermi 
level. Consequently, gating effect on the transmission, 
Fig. ||(b), and current. Fig. ^(b), is smaller compared to 



the top-site situation. 

For the hollow-site configuration there is no trace of the 
HOMO-LUMO gap in the I-V curve, and the I-V curve 
is almost ohmic in the wide range of voltages V < 2V 
(Fig. This is indicative of the metalization of the 
chemically bonded molecule. This should have general 
implications, the simplest being an obvious difficulty in 
gating such molecules. This relates well to the large 
observed values of conductance (close to a conductance 
quantum) through a sing le molecule, and its small tem- 
perature dependence |23|| . 

Finally, it is important to mention that in the present 
as well as other calculations, the LUMO is the closest 
molecular orbital to the Fermi level Ep and the maxi- 
mum gating effect is naturally expected when its energy 
is pulled down closer to Ep. This takes place at positive 
gating voltage, and not the negative one, as observed ex- 
perimentally . This is an apparent contradiction which 
needs to be resolved. 

C. Effect of contact geometry 

The tilting angle has a large effect on the LV curves 
of BDT molecules, see Figs. ^, The behavior of the 
BDT on the top site and on the hollow site is again rather 
different. The I-V curve for the hollow site remains ohmic 
for tilting angles up to 75° with moderate changes of 
conductance. Fig. ^(b). The variation of the current with 
the angle are much larger for the top site. Fig. ||(a). 
By changing from 5° to just 15°, one drives the I-V 
characteristic from one with a gap of about 2V to the 
ohmic one with a large relative change of conductance. 
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FIG. 5. Current-voltage characteristic and effects of gating 
and tilting with respect to the Au(lll) electrode surface on 
current through the BDT molecule on the hollow site: (a) 
effect of the gating, (b) effect of increasing tilt angle. Current 
is in units of Jo — 77.5fiA. 

Even changing d from 10° to 15° changes the conductance 
by about an order of magnitude. 

Finally, Fig. ^ illustrates the role of the electronic 
structure of the electrodes. The results described above 
have been obtained with only the s-states on Au atoms. 
We have also considered an sp- and spd-basis for Au. Al- 
though substantially increases the computing time, but 
the addition of p- and d-states brings about only mod- 
erate changes in current. Since the hybridization is dif- 
ferent for different cases, the current magnitude slightly 
varies for different basis sets. 

D. Possible origins of the observed gating effect 

With regards to the origin of the observed gating ef- 
fect, one can envisage that in dilute BDT-alkane solu- 
tion sandwiched between Au electrodes there may be 
two processes going on that significantly change conduc- 
tance. Since the BDT molecules are clamped by the ma- 
trix of alkane chains, they have to move with it. Note 
that the matrix of alkane chains is not in registry with 
the Au electrode. The BDT molecules have a nominal 
length of 7.2Aand are dissolved in the matrix of (CH2)5S 
alkanethiol insulating molecules with a nominal length of 
8.3A. Thus, the BDT molecule would appear as a dip on 
the surface of the matrix, and one or a few gold atoms 
can get into this dip during the deposition and bind to 
the end S. The geometry of this bond is uncertain, and 
the bond may well be stretched. In this case even a slight 
perturbation exerted on the SAM might lead to reconfig- 
uration of the bond resulting in large changes of conduc- 
tance. It seems reasonable to assume that the "domain 
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FIG. 6. Effect of electrode structure of Au(lll) electrodes 
on current through the BDT molecule on the top site. There 
is moderate difference between calculations using s-, sp-, and 
full spd-bases. Current is in units of Jo = 77.5^^4. 

walls" separating different patches in an alkane matrix 
move rather freely in the system, since it does not require 
much energy. The BDT molecules will follow the matrix 
and can either snap from a hollow site to a top site and 
back and/or change the tilt angle. Both processes may be 
accompanied by large changes in the conductance. Con- 
formational changes of the clamped BDT molecule are 
rather restricted, and the motion of the "domain walls" 
may be quite repeatable. One may wonder what causes 
the domain walls to move. As a possible reason, we sug- 
gest the presence of positive metal ions inside the or- 
ganic film, as a small concentration of electrode ions in 
a SAM is rather inevitable. Indeed, Au"*" strongly inter- 
acts with CeHg in the gas phase and forms an Au+— CeHg 
complex with a binding energy 2.65 eV, whereas neutral 
Au forms a Van-der-Waals complex with the binding en- 
ergy 90 meV |2j]. Even more likely is a formation of 
those complexes with thiophenes, which carry an elec- 
tric dipole. It is likely that a similar charged complex 
can form with BDT molecules in a SAM with those BDT 
molecules that have lost or are in poor contact with the 
gold substrate. It is also possible that a charged complex 
can be formed between the charged metallic ions (Au"*" 
or other electrode metals) and alkane chains. A small 
field in the organic film will then produce a tangential 
force on the ions, and this may trigger the domain wall 
motion when the pinning is weak. Additionally, since the 
packing of the film is not ideal (an organic film is usually 
a rather disordered patchwork of "grains"), the Maxwell 
force acting on the top Au electrode at finite drain voltage 
and/or electron wind force may trigger the domain wall 
motion, which may also require a combination of these 
factors. The second possibility would be a build-up of 
the interface charge, but apparently the replacement of 
oxide did not change the results. As mentioned above, it 
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is also difficult to explain the jump in conductance at a 
certain value of the gate voltage. 

IV. SUMMARY 

We have presented an extensive analysis of the 
electronic states and transport through the benzene- 
dithiolate molecule, which is the simplest conjugated 
molecule that forms a SAM and exhibits large changes 
of conductance under gating It shows that the ef- 
fect of gating strongly depends on the geometry of the 
molecule-electrode contact, and is maximal for the less 
coordinated top-site position. It is related to the sharp- 
ness of the peak in transmission, which corresponds to 
the LUMO on BDT and is close to the Fermi level Ep of 
an Au electrode. By the same token, the current is more 
sensitive to the tilting angle of the molecule when it is 
positioned on the less coordinated top site. It is worth 
mentioning that the fact that the LUMO is closer to Ep 
suggests that a positive gating voltage Vg should produce 
the larger effect, not a negative gate voltage, as was ob- 
served in [Q. This discrepancy should be addressed in 
the future. Binding on the highly coordinated hollow 
site naturally leads to large hybridization of the molec- 
ular states with electrode states, which become smeared 
out. Consequently, BDT molecule becomes metallized, 
i.e. the I-V characteristic becomes practically ohmic. A 
very small effect of gating is predicted for this geometry. 
In any case, it is difficult to explain the large observed 
effect of gating on BDT molecules in a slot with a width 
of only 1 nm by small voltage applied to the gate 4-5 
nm away. One should assume that there is either (i) a 
build-up of interface charge in the immediate vicinity of 
the slot opening, which is very sensitive to the gate volt- 
age; or (ii) small inhomogeneous electrostatic forces and 
resulting stresses on the SAM result in the reconfigura- 
tion of the film and, consequently, of a clamped inside 
BDT molecule with respect to the gold contacts. Both 
of those mechanisms have problems of their own in ex- 
plaining the experimental observations, as discussed in 
the text. Changes in random charges in the film and/or 
chemical composition of the BDT molecules (e.g. loss of 
end sulfur) are possible but are unlikely to be reversible. 
One needs to characterize the films better and vary the 
gate oxide thickness and other parameters of the system 
in order to confirm one of those mechanisms or suggest 
some other effects controlling the gating in slot geome- 
tries like the one used in Refs. ||7|,pl|. 

We thank J.H. Schon for extensive helpful conver- 
sations and sharing his data and preprints, and R.S. 
Williams for useful discussions. The work has been partly 
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